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Executive  Summary 


Impedance  matching  is  a  canonical  problem  in  electrical  engineering — design 
a  matching  network  that  transfers  maximum  power  from  a  generator  to  a 
load.  Although  the  loads  in  this  report  are  restricted  to  compact  antennas, 
the  mathematical  exposition  applies  to  any  lumped  load. 

Finding  a  “good”  matching  network  is  a  difficult  numerical  optimiza¬ 
tion  problem.  Therefore,  techniques  that  compute  bounds  on  the  matching 
provide  valuable  information.  For  example,  the  Fano  bounds  discussed  in 
this  report  bound  the  matching  performance  as  a  function  of  the  frequency 
band.  Therefore,  the  circuit  design  can  see  the  matching  performance  against 
the  bandwidth  to  make  design  decisions — without  having  to  solve  countless 
matching  problems.  Moreover,  the  Fano  bounds  provide  an  excellent  bench¬ 
mark  to  assess  the  performance  of  actual  matching  circuits — getting  “close 
enough”  to  the  best  bounds  is  “good  enough.” 

The  Fano  bounds  are  computed  by  maximizing  the  matching  performance 
under  inequality  constraints.  The  inequality  constraints  are  determined  by 
analytic  expansions  of  the  load.  Because  symbolic  manipulators  compute 
these  expansions  and  the  numerical  packages  can  solve  the  inequality  con¬ 
straints,  the  Fano  bounds  may  be  amenable  to  a  “hands-free”  computation, 
ffowever,  the  sticking  point  for  the  Fano  bounds  is  the  requirement  that  the 
load  be  given  in  the  Darlington  representation.  Therefore,  this  Phase  I  ef¬ 
fort  shows  that  the  Fano  bounds  can  automate  and  defers  the  Darlington 
computation  to  Phase  II. 
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1  Antenna  Systems 


The  key  idea  for  antenna  design  is  that  both  the  antenna  and  its  matching 
circuit  are  the  objects  under  design — not  the  antenna  with  the  matching 
circuit  as  an  afterthought.  This  point  is  developed  in  Zhu  and  Qi  [30]: 

The  “raw”  impedance  characteristic  of  an  antenna  . . .  does  not 
provide  a  ready  index  to  what  might  be  described  as  the  “broad¬ 
band  potential”  of  the  antenna. 

The  notion  of  an  antenna  system,  is  found  in  Gustafsson  and  Nordebo  [13]: 

The  bandwidth  of  an  antenna  system  can,  in  general,  only  be 
determined  if  this  impedance  is  known  for  all  frequencies  in  the 
considered  frequency  range.  However,  even  if  the  impedance  is 
known,  the  bandwidth  depends  on  the  specified  threshold  level  of 
the  reflection  coefficient  and  the  use  of  matching  networks. 

This  notion  of  an  “antenna  system”  is  unavoidable  when  designing  integrated 
systems  such  as  transmit /receive  modules  and  multiple-input /multiple-output 
systems.  For  example,  a  wideband  microstrip  circulator  connected  to  an  an¬ 
tenna  requires  a  new  design  methodology  [29] : 

that  allows  for  ...  a  single  matching  network  that  interfaces  be¬ 
tween  the  circulator  and  the  antenna  . .  . 

In  more  detail,  matching  networks  trade  isolation  between  the  transmitter 
and  receiver  against  power  transfer  between  the  transmitter  and  the  antenna, 
and  the  antenna  and  receiver.  More  generally,  the  limited  bandwidth  of  the 
circulator  forces  the  designer  to  trade  power  transfer  and  isolation  to  get  more 
bandwidth.  Consequently,  antenna  “tweaking”  is  a  poor  design  approach  for 
these  integrated  systems. 

Basic  to  the  performance  of  any  antenna  system  are  the  matching  net¬ 
works.  Classically,  these  matching  networks  are  designed  to  optimize  power 
transfer.  This  classic  design  objective  is  expanding  to  encompass  multiple 
design  objectives — increasing  bandwidth  [16],  [2],  increasing  channel  capacity 
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[18],  reducing  bit  error  rate  [27],  [6].  In  any  optimization  problem,  knowing 
global  “best  bounds”  provides  a  benchmark  to  assess  designs,  terminate  the 
numerical  optimizers,  and  end  the  interminable  “tweaking”  by  the  engineers. 
The  best  matching  performance  is  currently  computed  by  the  following  tech¬ 
niques: 

•  Analytic  matching  theory  of  Fano  [9]  or  Youla  [28], 

•  H°°  methods  of  Helton  [14] , 

•  Real-frequency  technique  of  Carlin  and  Civalleri  [5] , 

•  Analytic  bounds  of  Wohlers  [26] , 

•  Impedance  disks  of  Sokolov,  Babak  and  Babak  [24], 

A  complete  report  would  compare  these  methods  on  real  antennas.  Phase 
I  starts  this  process  with  the  Fano  Method.  The  immediate  objective  is  to 
implement  Fano’s  integral  expansions  in  a  symbolic  calculator.  The  long- 
range  objective  is  to  automate  the  Fano  Method  to  compute  matching  bounds 
starting  from  numerical  antenna  measurements. 

Section  2  makes  explicit  the  assumptions  behind  the  key  inequality  ob¬ 
tained  by  Fano  in  his  1950  thesis  [9] : 

1  —  Gx  Y  |  sx  | 2  • 

This  inequality  shows  the  transducer  power  gain  Gt  of  any  matching  net¬ 
work  is  bounded  by  the  Fano  reflectance  sx  extracted  from  a  Darlington’s 
representation  of  the  load.  Consequently,  the  Fano  reflectance  from  the  fixed, 
load  limits  the  best  possible  gain  achievable  by  any  lossless  matching  network 
selected  by  the  circuit  designer. 

The  gain-bandwidth  bounds  are  computed  by  evaluating  integrals  of  the 
form: 

/OO 

log(| sx(ju)\~l)du  =  -Af  -v T^pm. 

Section  3  reviews  a  convenient  mathematical  technique  to  evaluate  these 
integrals.  Appendix  A  examines  this  expansion  to  clarify  the  treatment  of  the 
DC  term.  Comparison  is  made  with  the  symbolic  calculators  used  throughout 
this  report  for  automating  the  Fano  computations. 
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Section  4  combines  the  integral  techniques  with  the  Fano  inequality  to 
develop  Fano’s  gain-bandwidth  bound: 

,°g(i  _  Gm)  <  r iog(|w(iw)r,)ljWi 

where  G'max  denotes  the  maximum  gain  attainable  over  the  frequency  band 
hi  =  [aAnin,  14m].  Section  5  works  out  the  Fano  bounds  restricted  to  a  nor¬ 
malized  resistor.  Section  6  drops  the  normalization  constraint  so  that  the 
Fano  bounds  for  antenna  matching  can  be  computed. 

A  Fano  analysis  for  a  microstrip  patch  antenna  and  ideal  dipole  is  under¬ 
taken  in  Sections  7  and  8,  respectively.  The  careful  development  provides  a 
more  flexible  approach  to  the  Fano  bounds.  Indeed,  the  literature  typically 
obtains  the  Fano  bound  from  a  system  of  equality  constraints.  For  example, 
one  system  has  the  form 


71 K 
~~2 
7 tK 

- I 

6 


( 


C(7n 


,  ,3  _  3  \ 

^max  '-‘•'mill ) 


log(kvCM|  )du, 


V  log(|s.YCM|  )du, 


where  K  is  standard  notation  for  the  normalized  bound  on  the  Fano  load 
over  the  frequency  band: 


min{log(|s.Y(jcn)|  x)  :  u  G  hi}  =:  nK/2. 


The  Fano  bound  is  computed  by  solving  for  these  equality  constraints  for  K 
and  then  backing  out  the  transducer  power  gain.  As  long  as  the  two  integrals 
are  transverse  with  respect  to  K,  a  unique  solution  exists.  This  approach 
breaks  down  as  additional  constraints  are  derived.  Unfortunately,  “solutions” 
found  in  the  literature  are  computed  by  tossing  all  but  two  of  the  equality 
constraints  to  obtain  a  sub-optimal  bound.  However,  the  Fano  bound  is  really 
constrained  by  a  system  of  inequalities.  In  terms  of  the  preceding  example, 
the  Fano  bound  is  computed  by  maximizing  K  subject  to  the  inequalities: 


7 tK 
ttK 

ir( 


;3 

'max 


^min  ) 

< 

_  ,  .3  \ 

^min  ) 

< 

log(|sYCM|  )du, 


U  log(|sYCM|  )du. 
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An  optimization  problem  with  inequality  constraints  is  a  special  case  of  the 
Kuhn-Tucker  theory  [17,  Section  9.4],  Because  the  Kuhn-Tucker  theory  was 
not  formalized  until  11  years  after  Fano’s  1950  thesis,  it  appears  that  the 
special-case  solutions  obtained  from  the  system  of  equalities  have  been  con¬ 
flated  with  the  solutions  to  the  general  system  of  inequalities.  In  this  report, 
the  antenna  examples  compute  the  Fano  bounds  by  carefully  working  through 
the  inequality  constraints. 

Two  major  problems  impede  a  “hands-free”  implementation  of  the  Fano 
bounds  for  real-world  antenna-matching  problems: 

1.  Integral  expansions  must  be  automated. 

2.  Numerical  antenna  measurements  must  be  mapped  into  its  Darlington’s 
representation. 


Section  9  discusses  both  problems  in  detail  and  sets  out  the  approach  for 
Phase  II.  The  symbolic  calculations  used  in  this  report,  coupled  with  re¬ 
stricting  the  load  to  antennas,  demonstrate  that  the  integral  expansions  can 
be  automated — albeit  a  non-trivial  piece  of  programming. 

The  problem  mapping  of  a  finite  number  of  antenna  measurements  to  a 
Darlington  representation  is  more  difficult.  For  example,  the  Fano  bounds 
for  some  loads  are  based  on  the  expansion  of  the  reflectance  sx(p)  at  p  —  oo. 
Consequently,  both  the  presence  of  the  transmission  zero  at  p  =  oo  and  its 
order  must  be  extracted  from  the  numerical  measurements — and  a  priori 
knowledge  of  the  antenna’s  high-frequency  behavior.  Appendix  B  demon¬ 
strates  that  a  Darlington’s  representation  can  be  obtained  from  sampled 
antenna  measurements.  The  quality  of  this  representation  can  be  enhanced 
by  incorporating  a  priori  knowledge  of  the  antenna — especially  its  behavior 
at  high  and  low  frequencies.  Finally,  progress  made  in  mapping  the  antenna 
measurements  carry  over  to  other  impedance-matching  methods,  especially 
the  H°°  methods  that  also  rely  on  knowledge  of  the  load  over  all  frequencies. 

The  notation  is  the  standard  scattering  formalism  and  illustrated  in  Fig¬ 
ure  1.  The  voltage  and  current  vectors 


v  = 


V\ 

V2 


1  = 


h 

12 


4 


of  the  2-port  determine  the  incident  wave  [3,  Eq.  4.25a],  [4,  page  234] 


a  = 


a  i 
0-2 


?V2; 


and  the  reflected  wave  [3,  Eq.  4.25b],  [4,  page  234] 


b  = 


bi 

b2 


=  1/2v  -  Ri" i}, 


?V2  J 


with  respect  to  the  normalizing1  matrix 


(1) 

(2) 


Consequently,  the  scattering  matrix  S  maps  the  incident  wave  a  to  the  re¬ 
flected  wave  a 


b  = 

'  6i  ' 

Sll 

Sl2 

Oi 

b2 

S21 

S22 

02 

Ssl. 


Figure  1:  The  2-port  scattering  formalism. 

The  scattering  matrix  determines  the  impedance  matrix  Z  as 

Z  :=  Rfl1/2ZR-i/2  —  (I  +  S)(I  —  S)~\ 

1Two  accessible  books  on  the  scattering  parameters  are  Baher  [3]  and  Balabanian 
&  Bickart  [4].  Baher  omits  the  factor  of  1/2  but  carries  this  rescaling  onto  the  power 
definitions.  Most  other  books  use  the  power- wave  normalization  [8] :  a  =  Rq  '  {v+Zoi}/2, 
where  the  normalizing  matrix  Z0  =  R0  +  jXq  is  diagonal  with  diagonal  resistance  Rq  >  0 
and  reactance  X0. 
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provided  /  —  S'  is  invertible.  Here  Z  denotes  the  normalized  impedance 
matrix.  If  the  iV-port  admits  an  impedance  matrix  Z, 

S  =  {Z-I){Z  +  I)-\ 

Table  1  summarizes  the  notation. 


Table  1:  Table  of  notation. 


Symbol 

Description 

J 

V-l 

V 

Complex  frequency:  p  =  a  +  joj 

OJ 

Radial  frequency 

f 

Frequency  (Hertz):  uj  =  2i r/ 

a 

Neper  frequency 

K 

Real  part:  5R[p]  =  a 

Imaginary  part:  ^[p]  —  u 

R 

Real  numbers 

C 

Complex  numbers 

C+ 

Open  right  half  plane 

jR 

Imaginary  axis 

r0 

Normalizing  resistance 

SL 

Reflectance  of  the  load 

sg 

Reflectance  of  the  generator 

sx 

Fano  reflectance 

S 

Scattering  matrix 

zL 

Impedance  of  the  load 

Zx 

Fano  impedance 

Z 

Impedance  matrix 
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2  Fano’s  Approach  to  Wideband  Matching 


Most  impedance  matching  techniques  look  into  the  load.  In  contrast,  Fano 
looks  out  from  the  load  to  compute  fundamental  limitations  on  wideband 
matching.  Specifically,  Fano  uses  Darlington’s  theorem  to  compute  wide¬ 
band  matching  bounds  over  the  class  of  all  lossless,  lumped  2-ports.  This 
section  reviews  the  lossless  2-port  theory  to  implement  Fano’s  view.  Figure  2 
shows  the  incident  waves,  reflected  waves,  and  power  flows  of  a  2-port  with 
scattering  matrix 

ip  =  + 

that  connects  a  generator  to  a  load. 


S(p)  = 


snip)  snip) 
S2lip)  S22ip) 


Figure  2:  Reflectances  and  powers. 

The  2-port  is  called  lossless,  provided  S  is  unitary  over  all  frequencies: 

S(ju)HS(jw)  = 

The  reflectance  looking  into  Port  1  is 

s i  =  —  =  S'i i  +  <Si2<Sl(1  —  s22Sl)  1<S2i  =:  T i (S,  Sl).  (3) 

CL\ 

Similarly,  the  reflectance  looking  into  Port  2  is 
&2  1 

=  —  =  S22  +  S2lS(?(l  -  SuSG)  s12  =:  ^(S,  SG).  (4) 

a  2 


1  0 
0  1 
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The  power  flows  are  determined  from  the  incident  and  reflected  waves: 


•  The  average  power  delivered  to  Port  1: 

Pi-=  5(M2-N2)  =  ^(i-|siR- 

•  The  average  power  delivered  to  Port  2: 

i(M2-N2)  =  -/V 

•  The  average  power  delivered  to  the  load  [11,  Eq.  2.6.6]: 

Pl-=  l(K|2-|6i|2)  =  M(l-|Sl|2). 

•  The  average  power  delivered  by  the  generator: 

Pt  kcl2  i-N2 

D  |  -1  19’ 

2  1 1  sqS\  p 

The  maximum  power  available  from  a  generator  is  defined  as  the  average 
power  delivered  by  the  generator  to  a  conjugately  matched  load  [11,  Eq. 
2.6.7]: 

-PG,max  :=  Pg\Si=s£  =  ~~~  (1  —  Isg|2) 

The  source  mismatch  factor  [11,  Eq.  2.7.17]  is 

Pg  (i  -  kcl2)(i  -  ki|2)_ 

Pg,  max  |1  —  SG-Sl|J 

The  maximum  power  available  from  the  2-port  is  defined  as  the  average  power 
delivered  from  the  network  to  a  conjugately  matched  load  [11,  Eq.  2.6.19]: 

p  . _  pi  _  . _  I  ^2 1  Sl=S2  |  /'1  _  I  |2\ 

L, max  •  *  L\sl=S2  •  2  |^2|  )• 

The  load  mismatch  factor  [11,  Eq.  2.7.25]  is 

PL  (l-|a£|2)(l-|52|2) 

Pl,  max  |1  SlS  2^ 
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The  basic  gains  are  itemized  [11,  page  213]: 


•  Transducer  power  gain: 


Gt  := 


Pl 


power  delivered  to  the  load 
maximum  power  available  from  the  generator 


•  Power  gain: 


Pl  power  delivered  to  the  load 
P\  power  delivered  to  the  network 


•  Available  power  gain: 

Pt,  may  maximum  power  available  from  the  network 

Gr a  •—  -  —  - 7 - 7 - . 

Pg, max  maximum  power  available  from  the  generator 


When  a  general  2-port  is  used  in  Figure  2 — be  it  lossless,  passive,  or  active — 
the  transducer  power  gain  Gt  is  provided  by  Pozar  [20,  Eq.  11.13]: 


Gt(sg,  S,  $l) 


1  So,  I2  - 

1  - 

|sg!: 

2 

1  - 

|  SL\ 

2 

lS2l| 

1  - 

"|2  | 

1  -  -S22Sl| 

(5) 


The  smallest  transducer  power  gain  over  a  frequency  band  11  =  [cumin ,  a>max] 
is  denoted: 


||GT(sG,*S,,SL)||f2-oo  :=  sup {GT(sG(jw),S(ju),sL(ju}))  :  w  e  O}. 

The  goal  of  impedance  matching  is  to  maximize  the  smallest  transducer 
power  gain: 


Impedance  Matching:  Maximize  the  transducer 
power  gain  over  a  given  a  class  S  of  lossless  2-ports: 

sup(II^t(sg,  S,  sl)||o,-oo  :  S  €  5}. 
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Impedance  matching  admits  several  equivalent  formulations.  Define  the 
power  mismatch  between  reflectances  Si  and  S2  as  [1]: 


AP(s1,s2) 


si  -  S2 
1  -  sis2 


When  the  2-port  is  lossless, 


Gt(sq,  S,  sl )  =  1  —  A P(sc,  si)2 


1  - 


sg  ~  s  1 


2 


1  -  SGS 1 


(6) 


Consequently,  the  power  mismatch  reveals  the  conjugate  matching  embedded 
in  the  transducer  power  gain.  Moreover,  a  lossless  2-port  admits  matching 
at  either  Port  1  or  at  Port  2. 


Lemma  1  [1]  If  the  2-port  of  Figure  2  is  lossless  and  |sisG|,  |s2Sl|  <  1, 
CrT(sG ,  sL)  =  1  -  A P(sG,  si)2  =  1  -  AP(s2,  sl)2. 

The  largest  power  mismatch  over  a  frequency  band  =  [o;min,  o;max]  is  de¬ 
noted: 


II  AP(sg,  Si)  ||  12,00  :=  sup  {AP(sG(jv),  Si(ju))  :u60}. 

Because  the  smallest  transducer  power  gain  is  the  largest  power  mismatch, 

||Gt(sg, -S', sL)\\n ,-00  =  1  -  ||AP(sG, si)||nj0O. 

Consequently,  maximizing  the  smallest  transducer  power  gain  is  equivalent  to 
minimizing  the  largest  power  mismatch.  Fano  [9]  uses  Darlington’s  theorem 
and  Lemma  1  to  untangle  the  load  from  the  matching  network. 

Theorem  1  (Darlington)  [26],  [3]  Any  rational  bounded-real  function  is 
realized  as  the  reflectance  of  a  rational,  lossless,  reciprocal  2-port  terminated 
in  a  1-ohm  resistor. 
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Darlington’s  theorem  provides  a  representation  of  the  load’s  reflectance 
sl(p)  as 

sL(p)  =  F i(SL,  0;  p), 

where  Sl(p)  is  a  2x2  lossless  scattering  matrix  terminated  in  a  1-ohm  resistor. 
Figure  3  shows  the  matching  2-port  connected  to  the  load  s i  (p)  given  by 
Darlington’s  representation.  It  will  be  convenient  to  call  the  reflectance  Sx{p) 
looking  in  from  the  1-ohm  resistor  the  Fano  reflectance.  The  transducer 
power  gain  Gt.x  from  the  generator  to  the  1-ohm  resistor  is 


Gt.x  — 


P 


x 


Pt  P 


x 


P 


G,max 


P< 


G,max  Pl 


—  GtGxi 


where  Gx  is  the  power  gain  from  the  input  of  the  load  to  the  1-ohm  resistor. 


IQ 


Figure  3:  Matching  2-port  and  Darlington’s  representation  of  the  load. 
Because  G\  never  exceeds  1, 


Gt.x  >  Gt-  (7) 

Because  both  2-ports  are  lossless,  Lemma  1  links  the  Fano  reflectance  sxip) 
to  the  transducer  power  gain: 

Gt.x  =  1  —  AP(sx,  0)2  =  1  —  |sx|2  >  Gt, 

where  the  zero  in  the  power  mismatch  implies  the  normalizing  resistor  is  also 
1-ohm.  Consequently,  the  system  loss  is  bounded  by  the  Fano  reflectance: 

(8) 


1  —  Gt  >  |-s.y|2 
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Because  the  matching  network  appears  only  in  the  gain,  the  Fano  reflectance 
bounds  the  performance  of  any  lossless  matching  network: 

1  -  IMIL  >  sup{||GT(sG,5,sL)||00  :  S  G  5}. 

This  bound  is  a  simple  application  of  Equation  8.  Fano  uses  the  analyticity 
and  structure  of  sx(p)  to  compute  the  more  sophisticated  gain-bandwidth 
bounds.  The  first-order  bounds  are  obtained  from  the  Bode  theorems. 
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3  Bode  theorems 


If  h(p)  is  a  rational  function,  the  relative  degree  of  hfp)  is  denoted 

degr[/i]  =  #  of  zeros  at  infinity. 

Equivalently,  if  the  relative  degree  of  h(p)  is  positive,  h(p)  admits  the  expan¬ 
sion  at  infinity  as 

%)  =  — +  ^r  +  ---  (p->  oo), 

p  pz 

where  the  first  non-zero  term  h-d  has  index  d  =  degr  [h] . 

Theorem  2  (Bode  Integral  [12])  Let  h(p)  be  a  real  rational  function  that 
admits  an  expansion  at  infinity  as 

Kp)  =  —  +  ^  (p->oo). 

p  pz 

If  g(p)  :=  (1  +  h{p))~l  has  neither  zeros  nor  poles  in  the  closed  right  half 
plane  C+  , 

7 T 

\og(\g(ju)\)du  =  - -h_i. 


Theorem  3  (Modified  Bode  Integral  [12])  Let  h(p)  be  a  real  rational 
function  that  admits  an  expansion  at  infinity  as 

h(p)  =  —  +  ^  +  ...  (p->oo). 

p  pz 

If  g(p)  :=  (1  +  h{p))~l  is  analytic  on  C+  and  has  zeros  pi,  pm  in  the 
open  right  half  plane  C+  , 


\og(\g(jv)\)duj 


7T  M 

--h-1  +  7T  y^/Pm, 

m=  1 
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Example  1  (Single  Pole)  Let  a6l  and  define 


1  1  a  a2  a3 

Kp)  = - -  =  -  +  ^  +  ^  +  ^r  +  --  -  (p  ->  oo) 


p  —  a  p  p2  p 3  pA 


The  requirement  that  the  integrand 


9a  ip)  = 


p  —  a 


1  +  h{p)  p  —  (a  —  1) 
be  analytic  on  the  closed  right  half  plane  forces  a  <  1.  Theorem  3  computes 


7T 


log(|  ga(jw)\)dw  =  --  +  7r  a. 


Figure  4  plots  the  integrand  when  a  =  1/10. 


o 

-o.i 

-0.2 

-0.3 

-0.4 

-0.5 

-0.6 

-0.7 

-0.8 

-0.9 

-1 


log( |g(/o> )|);  g(p)  =  l/(  1  +h(p)  y,h(p)=\!(p-  1/10) 


r 

Jo 


log(\g(juj)\)duj  =  -2/5  7r 


10 

co 


12 


14 


16 


18 


20 


Figure  4:  Bode  integrand  for  Example  1. 
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The  figure  also  shows  the  symbolic  evaluation  of  the  integral  computed 
using  the  code  listed  Table  2. 

Table  2:  MATLAB®  code  evaluating  Example  1. 


syms  p 

syms  w  positive 
h  =  1/ (p  -  1/10); 
g  =  l/(  1  +  h  ) ; 
gW  =  subs(g,p, j*w) ; 

gBIT  =  int (  log (  abs(gW)  )  ,  0  ,  inf  ); 


The  first  application  of  the  Bode  integral  theorems  mix  the  impedance  and 
the  scattering  formalism  for  convenient  evaluation  of  the  Fano  reflectance. 


Corollary  1  Let  z(p )  is  a  real  rational  function  with  an  expansion  at  infinity 
as 


If  9±(p) 


z(p)  =  —  +  ^y  +  ...  (p->  oo). 

p  pz 

(1  ±  t(p))^1  are  both  analytic  on  C+, 


z(jv)  ~  1 

z(ju)  +  1 


duj  =  —nz-i. 


Proof:  The  Bode  integral  theorems  can  be  applied  to  the  sum  of  the  loga¬ 
rithms.  If  z{p)  =  u(p)/v(p)  where  u(p)  and  v(p)  are  real  polynomials  that 
are  relatively  prime, 


1  =  vjp) 

±  ^  1  ±  z(p)  v(p)  ±  u{p) 

Because  u(p)  and  v(p)  are  relatively  prime,  there  can  be  no  cancelation  be¬ 
tween  v (p)  and  u(p)  ±  v (p) .  Therefore, 

Zeros  [g±]  =  Zeros  [u] 

Poles [g±]  =  Zeros [u  ±  v] 
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Let  pi,  . .  . ,  Pm  denote  the  zeros  of  v (p)  in  the  open  right  half  plane.  The 
assumption  that  g±(p)  is  analytic  in  the  closed  right  half  plane  allows  the 
Bode  integral  theorems  to  be  applied: 


log 


z(juj)  -  1 


du 


z(juj)  +  1 

/»oo 

l°g(| z{ju)  -  1| )du  -  /  log(\z(ju)  +  1| )du 


\og(\z(ju)  ~  1\  )dw  +  /  \og(\z(ju)  +  1|  )du 


log(\9-{jv)\)dw  +  /  \og(\g+(ju)\)du 


i  o 


1 0 


M 


-  lim  p(-z(p))  +  7T  V  pri 

/  p— >oo  ‘ 

m= 1 

M 

--  lim  p(+z{p))  +  7T  V  pm 

/  p— >oo  ‘ 

m=  1 

=  — 7T^_! 


Ill 


Example  2  The  expansion  of  z(p)  at  infinity  is 

1 

p  +  3  p  p2 


zip)  = 


1  3  9 

o  I  o  *  *  ' 


P * 


C p 


and  produces  the  reflectance 


s{p)  = 


where 


9+ip)  = 


zip)  ~  1  =  9+ip) 
zip)  +  1  g-ipY 

1  p  +  3 


z  +  1  p  +  4 
is  analytic  on  C+  with  the  single  pole  at  p  =  — 4  and 

1  p  +  3 


9- ip)  = 


z  -  1 


p  +  2 


(X) 


is  analytic  on  C+  with  the  single  pole  at  p  =  —2.  Consequently,  Corollary  1 
applies  and  computes  the  integrand  as  —n. 
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Figure  5  plots  the  integrand  and  displays  the  symbolic  computation  ob¬ 
tained  from  the  code  listed  in  Table  3. 


log(|s(/'co)|) 


Figure  5:  Bode  integrand  of  the  reflectance  of  Example  2. 


Table  3:  MATLAB®  code  evaluating  Example  2. 


syms  p 

syms  w  positive 
z  =  l/(  p  +  3) ; 

sW  =  simplify (  subs(  (z-l)/(z+l)  ,  p  ,  j*w  )); 
gBIT  =  int(  log(  abs(sW)  )  ,  0  ,  inf  ) 
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A  more  flexible  computation  splits  the  integrand  and  accepts  separate 
right  and  left  half  plane  analyticity. 


Corollary  2  Let  z(p)  be  a  real  rational  function  with  an  expansion  at  infinity 
as 

7  i  7.  o 

(p  -*•  oo). 


/  N  Z-l  ,  *-2  . 

•nw  =  —  +  -5-  + 

p  pz 


If  g±ip)  (1  ±  z(p))  1  is  analytic  on  C±, 

I  z(Ju)  -  1 


log 


z(juj)  +  1 


du  =  7TZ—1  +  7T  ^  Pm, 


where  the  pm ’s  are  the  poles  of  z(p). 


Proof:  The  Bode  integral  theorems  apply  to  the  sum  of  the  logarithms.  Let 
z  —  u/v  be  real  polynomials  that  are  relatively  prime. 

1  v 

9z t  |  31  ' 

1  ± Z  V  ±u 

Because  u(p)  and  v(p)  are  relatively  prime,  there  can  be  no  cancelation  be¬ 
tween  v(p)  and  u(p)  ±  v(p).  Therefore, 

Zeros  [g±]  =  Zeros  [u] 

Poles[#±]  =  Zerosfw  ±  v] 


Let  pfi,  ...denote  the  zeros  of  v (p)  in  the  open  C±.  The  Bode  integral 
theorems  adapt  to  g~{p)  using  the  observation  that 


g*(p) 


9-(~p) 


1 

1  -  z(-p) 


is  analytic  on  the  open  right  half  plane.  The  Bode  integral  theorems  compute 


log\g-(ju)\)dw  = 


1 og\g-{-jv)\)duj 


~  \  -77  lim  p(-z(-p))  +  7rVp; 

I  Z  p—>oo  z — * 

-  p™} 

7T  _ 

+77 --1  -  n  z^pm- 
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Applied  to  the  reflectance,  the  Bode  integral  theorems  compute 


log 


z(ju)  -  1 


du 


z(juj)  +  1 

/»oo 

log(\z(ju)  -  l\)du  -  /  log(\z(ju)  +  l\)du 


/  log(\z(Ju)  -  1|  )duj+  \og(\z(ju)  +  1|  )dw 

'0  ^0 
/*  oo  /»oo 

/  log|^-C7'w)|)d£j+  /  log(|g+(jca)|)du; 

'o  ^0 

7T  \  ^  _ 

-  <!  +2^-i  -  n  Z^Pm 


Pr 


■^-1  +7T^p+} 

-7T^_i  +  VT 


III 


The  Bode  integral  theorems  apply  to  the  scattering  integrand  when  z(p) 
is  unbounded.  The  next  result  only  needs  one  pole  condition  and  admits 
several  useful  computations. 

Corollary  3  Suppose  z(p)  is  a  lumped  impedance — a  positive  real  rational 
function — that  admits  an  expansion  at  infinity  as 

z(p)  =  zip  +  z2p2  + ...  (p  — >  oo). 

Observe  that  the  reflectance 


s(p)  :  = 


zip)  ~  1 

z(p)  +  1 


is  analytic  on  C+.  Let  pi,  ...,  pm  in  C+  denote  the  zeros  of  sip)  in  C+. 
Then 


log(|st7a;)|)dcj  =  —  7T  lim 


P 


M 


oo  z{p)  ~  1 


+  7T  y^Pm- 


m=  1 
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Proof:  Write 


g(p)  ■= 


zjp)  ~  1 

z{p)  +  1 


1 

1  +  h{p) 


2 

zip)  -  1 ' 


Observe  that  if  h(p)  =  u (p)  jv (p) , 


u(p)  —  v(p) 
u(p )  +  vip) 


Because  u(p)  and  v(p)  are  relatively  prime,  (u(p)  +  v(p))  and  (u(p)  —  v(p)) 
are  relatively  prime.  Therefore,  there  are  no  pole-zero  cancelations  between 
this  numerator  and  denominator  of  g(p)  so  that 


J  Zeros [g\  =  Zeros [u  —  v] 
1  Poles  [g]  =  Zeros  [u  +  v] 


If  g(p)  is  analytic  in  the  closed  right  half  plane,  Theorem  3  computes 


logi\g(ju)\)du 


M 

-  2^-1  +  71  y^/Pm, 

m=  1 


Observe  that 

/i_ i  =  lim  ph{p) 

p — >oo 

to  get  the  final  form.  Ill 


2  lim 

p — >oo 


V 

zip)  -  1 


Example  3  If  z(p)  =  p, 

/  \  _  zip)  -  1  _  p  -  1 

S  P  zip)  +  1  P  +1 


has  a  pole  at  —  1  and  a  zero  at  p\  =1.  Because  \s(juj)  =  l,  the  integral  of 
l°g(|s(j£j)|_1)  is  zero.  Corollary  3  computes 


log(|s(jO)|)dw  = 


-7r  lim 

p — >OC 


P 


z{p)  ~  1 


+  7rpi  = 


-7r  lim 

p—>oo 


P 


p  —  1 


+  7T  =  0. 
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Example  4  If  zip)  =  p  +  p2 


s(v)  =  zM  ~ 1  =  p2+p~^ 

zip)  +  1  p2  +  p  +  1 

has  poles  in  the  left  half  plane  and  zeros  in  the  right  half  plane: 


Poles 

Zeros 

-0.5000  -  0.86601 

0.6180 

-0.5000  +  0.86601 

-1.6180 

A  symbolic  computation  of  the  integral  produces 


log(\s(ju)\)du  =  7T 


2 

V5  +  1 


=  7 r  x  0.6180. 


Corollary  3  computes 


\og(\s(ju)\)dw  = 


-7T  lim  — ^ - 

p^oo  x:(p)  —  1 


M 

+  7T  ^  pm 
m= 1 


=  0  +  7r  x  0.6180. 


Table  4:  MATLAB®  code  evaluating  Example  4. 


syms  p 

syms  w  positive 

z  =  p  +  p~2 ; 

s  =  simplify( (z-l)/(z+l) ) ; 

sW  =  simplify(  subs(  s  ,  p  ,  j*w  )); 

sBIT  =  int(  log(abs(sW))  ,  0,  inf  ); 
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4  Gain-Bandwidth  Bounds 


The  Bode  integral  theorems  can  be  coupled  with  the  Fano  bound  of  Equa¬ 
tion  8  to  produce  the  classic  gain-bandwidth  bounds: 

1  —  Gt  >  |  sx  | 2 

log(l  -  Gt)  >  21og(|sx|) 

-log(l-Gr)  <  -21og(|sx|) 

<=>  log(!  -  gt)  <  log( | sx re¬ 
integrating  this  inequality  produces  the  general  gain-bandwidth  integral  con¬ 
straints: 


1 

2 


log(l  -GT{ju))duj<  log(|s.Y(jn;)|  )du. 


As  a  special  case,  the  classic  gain  profile 


GT(juj) 


df max  <  UJ  <  U>2 

0  otherwise 


simplifies  the  inequality  to 

~U2  „  ^  l0g(l  ~  dr  max )  <  [  log(  |  {ju)  |  _1  )du. 
^  J  0 


(9) 


The  Bode-Fano  theorems  bound  the  integral  using  the  parameters  of  the 
Fano  load  sx(p)-  Consequently,  the  load  itself  bounds  the  maximum  gain. 
For  example,  suppose  the  load  is  an  inductor  in  series  with  a  resistor  as  shown 
in  Figure  6.  The  impedance  looking  in  from  the  resistor 


zx(p)  =  pL\  +z2(p) 


has  reflectance  normalized  to  the  resistor 


sx(p) 


zxip)  -  R 

zxip)  +  R 
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If  z<2  (p)  is  bounded  as  p  — >  oo,  Corollary  3  computes 


log(|5-Y(i^)|  )du=-K 


R 

U 


M 

~  n  ^  ^  Pm  j 

m=  1 


where  the  pm’ s  are  the  zeros  of  -Sx(p)  in  the  open  right  half  plane.  Combined 
with  the  classic  gain  of  Equation  9,  a  Bode  gain-bandwidth  bound  is  obtained 
[20,  page  262]: 

--7T  iogt1  -  hr  max  )  <  /  log(  |  SX  (jw)  \  ^  )du  <  ^ . 
z  Jo 

With  Acu  =  27tA /,  the  gain-bandwidth  bound  is 

—A/  log(l  —  Gmax)  <  y~ 


and  is  plotted  in  Figure  7. 

Now  consider  the  load  of  Figure  8.  The  impedance  looking  in  from  the 
resistor  [9,  pages  71-73] 

zx(p)  =  pL\  +  ( pC2  +  V2(p))~1 

has  reflectance  sx(p )  =  (zx(p)  —  R)(zx(p)  +  -R)-1  normalized  to  the  resistor. 
If  y^p )  is  bounded  as  p  — >  oo,  Corollary  3  computes 


log(|s.Y(jA)|  )du 


+■ 


nR 

17 


■  7T 


M 

'y  ^  Pm  ■ 

m=l 


23 


-A/log(l-r;  mJ=  RL;  (ohms/|.i  H) 


Figure  7:  Gain-bandwidth  for  the  RL  load. 


Although  there  is  an  additional  capacitor,  this  bound  is  exactly  the  preceding 
RL  bound.  Fano’s  inspired  solution  was  an  adroit  selection  of  analytic  win¬ 
dows  to  exploit  the  hidden  structure  of  the  load.  In  this  case,  Fano  observed 
that  Darlington’s  representation  of  the  load  forces  two  transmission  zeros  at 
p  =  oo.  These  transmission  zeros  force  the  expansion  at  infinity  to  have  the 
form 


log (sx{p)  X) 


/too  /too  /too  /too 

All  Alo  jT\a  Air 

—  +  -i-  +  ^-  +  Hb  +  .. 


p 


p° 


jr 


p° 


> 


where  the  first  two  terms  depend  only  on  the  load.  The  Bode  integral  really 
extracts  Af°: 


log(|s.v(jA)|  1)dw  =  +‘^Af> 


M 

^  ^  ^  Pm  ■ 

m= 1 
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L\ 


Figure  8:  Low-pass  RLC  load. 


Fano  uses  the  quadratic  window  to  extract  /hf: 

r°°  7t  ^ 

/  w2log(|  sx{ju)\-l)du  =  -~Af  +  n^pll. 

m=  1 

More  generally,  if  the  reflectance  has  an  N  transmission  zeros  at  infinity,  the 
expansion  at  infinity  has  the  form  [9,  Eq.  9] 


/too  /too  400  400 

log(Mp)-1)  =  IP  +  ~  +  ~pt  +  '  ' '  +  ^0TT  + 


where  the  leading  odd  coefficients  depend  only  on  the  lossless  network  of  the 
load.  The  corresponding  windowed  integrals  are 


^logd^wr1)^  =  (-i)B-  a~+1  - 


2  n  +  1 


J2Pr- 


for  n  —  0, 1, . .  . ,  At.  The  next  section  works  out  the  gain-bandwidth  bounds 
for  the  RLC  load  in  detail. 
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5  RLC  Fano  Bounds  (Unit  Resistor) 


Fano’s  gain-bandwidth  tradeoffs  for  the  low-pass  RLC  load  are  derived  in 
detail  [9,  pages  70-73].  Figure  9  shows  the  Darlington  representation  of  the 
load  and  the  Fano  reflectance  looking  in  from  the  1-ohm  resistor. 


Figure  9:  Low-pass  RLC  load  with  1-ohm  normalization. 


The  lossless  part  of  the  RLC  load  is  a  low-pass  ladder  of  degree  2  with 
scattering  matrix  [9,  Fig.  5]: 


Sxip)  = 


p(pL\  C‘2  +  Li  —  C'2)  2 
2  -p(pLiC2  +  C2  -  L\  ) 

L1C2  +  p(Li  +  C2)  +  2 


Consequently,  there  are  two  transmission  zeros  at  p  =  00.  The  impedance 
looking  in  from  the  resistor  [9,  pages  71-73] 

zxip)  =  pLi  +  ( pC2  +  V2(p))~1 


produces  the  Fano  reflectance 


sx  (p)  = 


Zxip)  -  1 


ZX(p)  +  1 

normalized  to  1-ohm.  The  expansion  of  the  integrand  at  infinity  is 


400  400  4°o  400 

logMp)"1)  =  ^  +  ^-  +  ^  +  ^  + 

p  p6  p4 

2  1  f  2  2  1 


—  +  —  (— - —  \  +  —2y +  o\p~ 5] 

pin  +  p3  \  3 L\  L\C2  I  +  p4  L\Cl  +  U[P  J’ 
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provided  the  admittance  y2  (p)  is  bounded  as  p  — >  oo.  The  expansion  shows 
that  the  load’s  transmission  zeros  push  the  matching  network  into  the  higher 
order  terms.  The  two  transmission  zeros  at  infinity  give  the  Fano  equations 
[9,  Table  1]: 

poo  _ _ 

J  log(| Sx(j^)\^)d(jJ  =  -Af  -TT^Pm 

poo 

j  ^2log(|  sx{ju)\~x)du  =  ~-A™  +  -J2p3m, 

where  the  pm' s  are  the  zeros  of  Sxip)  in  the  open  right  half  plane.  Denote 
the  low-pass  frequency  band  as  fl  :=  [0,  cuc].  Bound  the  performance  of  the 
reflectance  sx(p)  on  Q  as 

min{log(|s.Y(jca)|“1)  :  uj  G  fl}  =:  ttK/2. 

Consequently,  the  reflectance  sx(p)  is  bounded  on  the  entire  frequency  axis 
as 

i/i  /•  m-i\  ^  f  nK/2  \lu\  <  luc 
MkvMi [)  >  1 0 

Inserting  this  lower  bound  into  the  two  integrals  produces  the  inequalities 

poo 

uciiK/2  <  /  log(|s.Y(ja;)|“1)(icc; 

Jo 

poo 

lv^ttK/6  <  /  uj2log(\sx(juj)\~1)duj. 

■Jo 

Combining  these  inequalities  with  the  integral  equations  gives  [9,  Eq.  25,  26] 

uj<  <  Af  -  2  ^2  Pm 

y \K  <  -3A?  +  2J2pL 

The  problem  is  to  maximize  K  subject  to  these  Fano  equations.  It  was  Fano’s 
insight  that  equality  is  simultaneously  achieved  at  maximum: 

ucK  =  Ap  -2y/Pm 

<4K  =  -3AP  +  2^pI- 
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Because  the  complex  zeros  pm  come  in  conjugate  pairs,  only  the  real  part  of 
the  equations  need  to  be  retained: 


ucK  =  7L~-2^\m 

UJ3cK  =  -3A£  +  2  -  3UmUm, 


where  complex  frequency  p  =  u  +  jv  determines  that  $ft[p3]  =  u3  —  3 uv2. 

To  formulate  maximizing  K  as  a  standard  optimization  problem,  max¬ 
imize  the  first  equation  while  enforcing  equality  between  both  equations. 
That  is,  maximize 

^400  2 

K  =  h( u,  v)  =  — - 'V'  um 

Uc  UJC  ^ 

subject  to  the  constraints  u  >  0  and 


0(u,v)  = 


A? 

Uc 


+  3 


cu3 


^  Um  o  ^  3umVm  0* 


U, 


Ur 


By  the  Lagrange  multiplier  theorem,  if  h( u,  v)  has  a  local  extremum  under 
that  constraint  <7(11,  v)  =  0,  the  gradient  of  h  is  normal  to  the  tangent  space 
of  the  contours  of  the  constraint  function  g  [17,  Section  9.8]: 


V/'t  =  XV  g. 


In  this  case, 


Vh 


_6_  T  u2 

a,3  [0 


+ 


V 

2uv 


XV  g. 


Because  u  >  0,  the  v  component  of  the  Lagrange  equation  forces  v  =  0. 
Then  the  u  component  of  the  Lagrange  equation  simplifies: 


That  is,  all  the  iTs  are  equal.  Equivalently,  at  an  optimum  value  of  K,  there 
is  a  single  positive  zero  u0  that  must  satisfy  [9,  Eq.  27,  28] 


g(u) 


0. 
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The  single,  real  root  u0  of  g(u)  =  0  is 

u0  =  Q-lDl/*  -2uj2cD~1/\ 

where 

D  =  54  A^uj2c  +  162A^  +  6^48^  +  81  AfV  +  486 Afuj2cAf  +  72V  Af' 2 . 
An  optimal  value  of  K  is  linked  to  the  bandwidth  lvc  as 

K°»  = - ^^75 - '  (10) 

Taking  the  maximum  of  Equation  8  gives 

1  -  ||Gr||n,-oo  >  Ilsx||n,oo  =  exp(-vr Kopt). 

That  is,  the  transducer  power  gain  is  bounded  from  above  as 

1  -  exp(-vr/lopt)  >  ||GT||n,-oo-  (H) 

Figure  10  plots  this  upper  bound  on  the  transducer  power  gain  as  a  function 
of  bandwidth  for  nominal  values  of  the  inductor  and  capacitor.  The  plot 
shows  that  no  matching  circuit  can  deliver  more  than  half  power  to  the  load 
(Gt  =  —3  dB)  over  the  frequency  band  0-0.32  MHz. 

To  plot  the  Voltage  Standing  Wave  Ratio  (VSWR)  against  bandwidth, 
observe  that  the  input  reflectance  si  links  to  sx  as 

kiCMI2  =  1  -  GT(ju)  >  \sX(ju)\2 

because  the  1-ohm  normalization  sets  Sg  =  0.  Taking  the  maximum  yields 


INkoo  =  sup{|si(jo;)|  :w  6  fi}  >  ||s.v|koo  =  exp(-7rA'/2). 
Because  the  VSWR  is  an  increasing  function, 

VSWR(|k|koo)  >  VSWR(||sx||nj0o)  =  j  +  CXP< ~7r^v/^ 

1  —  exp(— irK/2) 

Consequently,  the  optimum  K  value  provides  the  lower  bound 


VSWR(||31|k,oo) 


1  +exp(-7r/lopt/2) 

1  -  exp(— 7rA'opt/2)  ’ 


(12) 


Figure  11  plots  this  lower  bound  as  a  function  of  bandwidth.  The  plot  shows 
that  no  matching  circuit  can  be  constructed  that  has  a  VSWR  less  than  3 
over  the  frequency  band  0-0.24  MHz. 
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Gain-Bandwidth  Bounds  for  low-pass  RLC 


A/(MHz):I  =50  pH;  <^=0.01  pF;i?=l  ohms 

Figure  10:  Upper  bounds  on  gain  as  a  function  of  bandwidth  ooc  =  27rA/. 


A/(MHz):i=50  nH;Cj=0.01  nF;/?=l  ohms 

Figure  11:  Lower  bounds  on  VSWR  as  a  function  of  bandwidth  c vc  =  2nAf. 
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6  RLC  Fano  Bounds  (Arbitrary  Resistor) 


Up  to  now,  the  load’s  resistor  has  been  normalized  out  of  the  analysis.  For 
practical  systems,  the  normalization  is  typically  set  to  50  ohms  while  the 
resistive  part  of  the  load  is  determined  by  the  physics  of  the  system.  Figure  12 
shows  a  matching  network  with  the  generator’s  side  normalized  to  i?0  =  50 
ohms  while  the  load  has  R  =  300  ohms  hanging  on  the  end. 


Figure  12:  Non-normalized  low-pass  RLC  load. 


Fano  [9]  normalizes  by  a  transformer  as  shown  in  Figure  13.  Looking 
into  the  transformer  from  the  matching  network,  the  transformer  has  a  chain 
matrix  [3,  page  20] 


T 


n  1  0 
0  n  \  ’ 


with  winding  number  n  =  \J Rq/R  to  transform  the  Rq  into  the  load’s  resis¬ 
tance  R  =  Ro/n2. 


Figure  13:  Low-pass  RLC  load  normalized  to  50  ohms. 
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Let  Zxip )  denote  the  impedance  looking  toward  the  generator  from  the 
input  port  of  the  transformer: 

zxip)  =  pLi  +  ( pC2  +  y2(p))_1. 

Looking  into  the  transformer  from  the  50-ohm  load  gives  the  impedance 

zxip)  =n2zxip )  =  J[Zx(l o). 

Normalizing  to  Rq,  the  reflectance  is 


zxip)  -  Ro  _  zx ip) / R  ~  1 

Zx  ip)  +  Ro  Zx  ip)/R+  1  ' 


(13) 


There  are  two  transmission  zeros  at  p  —  oo.  The  expansion  of  the  integrand 
at  infinity  exhibits  this  fact: 


log(sjr(p)  x) 


p  p 3  pA  p5 

2R  ,  1  f  2H?  2 R  |  ,  1  2 RyAp)  ,  _5 

Up  V\3L?  L\CP  p‘  L\(%  +  9 


provided  the  admittance  y2  ip)  is  bounded  as  p  — >  oo.  The  leading  coefficients 
now  contain  the  load’s  resistor.  Equation  10  computes  the  optimal  K  value 
assuming  normalization  by  Ro  instead  of  1  ohm.  Consequently,  Equations  11 
and  12  compute  the  upper  bound  on  the  transducer  power  gain  and  lower 
bound  on  the  VSWR — all  normalized  to  R0.  Figures  14  and  15  show  the 
corresponding  Fano  bounds  for  this  low-pass  RLC  circuit. 
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VSWR 


Gain-Bandwidth  Bounds  for  low-pass  RLC 


Figure  14:  Upper  bounds  on  gain  as  a  function  of  bandwidth  ujc  =  27rA /. 


A/(MHz):I  =50  Cj-0.01  nF;i?=300  ohms 

Figure  15:  Lower  bounds  on  VSWR  as  a  function  of  bandwidth  c vc  =  2nAf. 


33 


7  Microstrip  Patch  Antenna 

A  Fano  analysis  of  a  microstrip  patch  antenna  [10]  is  based  on  the  equivalent 
circuit  of  Figure  16.  For  reference,  Figures  17  and  18  plot  the  impedance 
and  reflectance  of  this  circuit  near  resonance. 

C4= 11.7  pF  Z3=15.6  nH 

o - )| - nrrr- 

C,=I4.07| 

0 - 


Figure  16:  Microstrip  patch  antenna — equivalent  circuit. 

The  Fano  analysis  starts  by  computing  the  transmission  zeros.  Let  Sx(p) 
denote  the  scattering  matrix  of  the  Darlington  representation  of  the  load. 
Port  1  looks  in  from  the  resistor.  Port  2  is  the  input  to  the  load.  The 
transmission  zeros  are  zeros  of  the  (1,2)  element  of  scattering  matrix  Sx(p)' 

2L1Cip2 

Sx,l2(p)  -  C2LlL3C4pi  +  ( 2L1C4  +  C2L1  +  L3c4)p 2  +  (C2L1C4  +  L3ilC4)p3  +  (C4  +  iii)p  +  1  ' 

This  rational  function  has  two  transmission  zeros  at  p  —  oo  and  two  trans¬ 
mission  zeros  at  p  —  0.  Rather  than  undertake  this  tedious  algebra,  a  more 
transparent  approach  extracts  the  transmission  zeros  directly  from  the  circuit 
diagram. 


Theorem  4  [25,  pages  165-167]  If  either  a  ladder’s  shunt  impedance  is  zero 
or  its  series  admittance  is  zero,  then  the  ladder  has  a  transmission  zero. 


Figure  16  shows  that  the  shunt  impedance  shorts  at  p  =  0  and  at  p  — 
oo.  Likewise,  the  series  admittance  is  an  open  circuit  at  both  p  =  0  and 
p  =  oo.  Theorem  4  identifies  these  four  cases  as  transmission  zeros.  Because 
the  degree  of  the  circuit  is  4,  these  are  the  only  transmission  zeros.  The 
reflectance  looking  in  from  the  resistor  is 

_  zx{p)  ~  Rq  _  zx{p)/R  -  1 
S A  P  zx (p)  +  Ro  zx (p)/R+  1  ’ 
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Micro  Strip  Patch  Antenna 


MHz 


Figure  17:  Microstrip  patch  antenna — impedance. 


using  the  notation  of  Equation  13  to  handle  the  resistor  R.  Expanding  at 

p  —  0: 


log  (sX(p)  X) 


A°0  +  A°lP  +  A°3p 3  +  A°y  . . . 

.  Li  2L2(3R2C2  +  3R2C±~L1)  3 

— J7T  +  2— P - — - p 


2z2(p)clL‘l 

R 


p4  +  0[p~% 


provided  z2{jp)  is  bounded  as  p  — >  0.  The  two  transmission  zeros  at  p  =  0 
push  the  load  z2 (p)  to  the  higher  order  terms.  These  two  transmission  zeros 
force  the  Fano  integral  constraints  [9,  Table  II] 


/•OO 

j  a;"2  log(|Sx  Wr1)^  =  -A\  -  tt  (14) 

/»oo 

J o  logd.xWr1)^  =  +  -  ^p“3,  (15) 
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Figure  18:  Microstrip  patch  antenna — reflectance. 


where  the  pm’s  are  the  zeros  of  Sx(p )  that  lie  in  the  open  right  half  plane. 
Expanding  at  p  =  oo: 


log (sx(p)  X) 


Aoo  Aoo  Aoo  Aoo 

.  ./J.  1  /I Q  -ii/1  -/Ac 

-4“  +  —  +  ^-  +  ^-  +  — ~  + 


-J7T  + 


p  p *  p° 

2  wl  2  3L1R2C2  +  3R2C2L3- 
C^R  Xp~  3  C2L1L3R3 


2^2 (p) 


P 


x 


1 

p3 


provided  z2  (p)  is  bounded  as  p  — >  00.  The  two  transmission  zeros  at  p  =  00 
push  the  load  z2(p)  to  the  higher  order  terms.  These  two  transmission  zeros 
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force  the  Fano  integral  constraints  [9,  Table  II] 

/OO 

l0g(| SxU^l^duJ  =  - A °°  -  TT^Pm, 

poo 

/  u i2  ]og(\sx(ju)\~1)du  =  —A?  +  ^rJ2p™- 


(16) 

(17) 


VT  7T 

-A2°  +  - 
2  3  3 

Again,  the  pm' s  are  the  zeros  of  sxip )  that  lie  in  the  open  right  half  plane. 

Denote  the  frequency  band  as  D  :=  [a;min,  cnmax].  Bound  the  performance 
of  the  reflectance  sx(p)  on  D  as 

min{log(|sx(jn;)|_1)  :  w  £  0}  =:  nK/2. 

Consequently,  the  reflectance  S\  (p)  is  bounded  as 


log(|s.Y(jw)|  )  > 


7Th/2  Ctlmin  <  |cu|  <  cumax 
0  otherwise 


Inserting  this  lower  bound  into  the  Fano  integrals  produces  the  following 
inequalities: 


7 TK 
~2 
nK 

- I 

6 

7 tK 

- I 

2 

7 tK 

- I 

6 


( ^max  ^min ) 


< 


■  , 3  _  .  ,3  \  < 
“'max  “-'min )  — 


Jo 


Amin  Anax) 


< 


,  3  _  3  \  / 

°“min  mnaxj  — 


log(|s.Y(jA)|  )du, 


v  log(|sYCM|  )duj, 


U  log(|sA-(jw)|  )du, 


V  '  log(|sx(jA)|  )du. 


Denote  the  zeros  of  sx(p )  as  pm  =  um  +  jvm.  Because  the  sums  are  real,  the 
Fano  integral  constraints  force  the  equalities 


3?  Pm 

^  l^2  'Pm 

3ft 


{^Pm 

K[Z^3]  =  J2 


^  ^  Um-i 

'y  ^  um  3 UmVmi 

E  JJm 

7/2  I  ? 

um  '  um 

Um{uh  -  3t>gJ 

(Um  +  VmY 
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Combining  these  inequalities  with  the  Fano  constraints  (Equations  14,  15, 
16,  17)  produces  the  following  inequalities: 


Is.  ((Uj, 


^min  ) 
,3 


K(u; 
K(u~3 


min  ^max ) 


<  A\ 

<  A\-2 


2  ^  ^  um, 


^(^max-^min)  <  +  2  Vum(u, 


\  ^ 


S?r  1 

m  ^  m)'i 


*  ^  0/2 


,.-3  ^ 

mm  ^max  ^ 


7/2  I  ^.2  ’ 

m  1  m 

<  -3^  +  2^^^l 


«+^m)3 

The  goal  is  to  maximize  K  subject  to  the  inequalities  determined  by  the 
zeros  in  the  right  half  plane.  Define  the  functions 


:=  UoV, 


LUr\ 


( ,  ,3  _  /  .3  \ 

l^max  ^min 


Mu,  v 
M  u,  v 

M  u,  v 
Mu,  v 

Maximizing  A'  is  equivalent  to  maximizing 


2  ^  ^ um  (  > 


1{-3A?  +  2j2umK-3v2m)}, 


(k'min  ^max) 


A? 


Ur 


( ,  ,  3  _  3  \ 

V^min  ^max 


‘|-3A§  +  2^ 


^  +  V‘n 

~  3<) 

(«m  +  Vm)3 


min{A0(u,v),  hi(u,v),  h2{ u,v),  h3(u,v)} 

subject  to  the  constraints  that  u  >  0.  Because  sxip )  is  a  real  rational 
function,  its  zeros  pm  in  the  right  half  plane  come  in  conjugate  pairs.  Coupled 
with  the  observation  that  the  hk(u,  v)’s  are  even  functions  in  v: 

M u,  -v)  =  hk( u,  v), 


the  optimization  domain  can  be  restricted  to  u  >  0  and  v  >  0.  The  following 
Lagrange  multiplier  is  a  classical  approach: 


max{/i0(u)} 

subject  to  the  constraints  that  u  >  0,  v  >  0  and 


0i(u)  :=  hi  (u)  ho  (u) 

02  (u)  :=  Mu)  -  Mu) 
03  (u)  :=  Mu)  -  Mu) 
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0, 

0, 

0. 


The  Lagrange  multipliers  take  the  form  [17] 


V/io  —  AiV<7i  +  A2  V  Q2  +  A3  V 


or  equivalently, 

(1  +  Ai  +  A2  +  A3)V/ro  =  AiV/ii  +  \2Vh2  +  A3V/13. 

Although  a  symbolic  manipulator  can  solve  this  algebraic  system,  and  the 
pTO’s  do  decouple,  the  multipliers  Ai,  A2,  A3  and  the  number  of  pm' s  are 
confounding  factors.  Fortunately,  the  long  and  difficult  algebraic  analysis 
is  evaded  by  exploiting  the  insights  gained  by  the  more  robust  approach  of 
simply  maximizing  the  function 

K  =  min{/r0(u, v),  /ii(u,v),  h2{ u,v),  h3( u,v)} 

subject  to  u  >  0,  v  >  0,  and  a  fixed  number  of  pm’s.  For  the  rest  of  the 
discussion,  it  will  be  more  convenient  to  use  the  same  inequalities  for  the  u’s 
and  the  v’s.  Accordingly,  this  approach  is  developed  and  the  insights  gained 
can  then  be  applied  to  the  Lagrange  multipliers. 


39 


Figure  19  shows  a  lower  bound  on  VSWR  as  a  function  of  fractional 
bandwidth  while  fixing  the  number  of  pm' s  to  1.  Specifically,  the  function 

K  =  mm{h0(ui,vi),  h^u^vi),  h2(u1,v1),  h3(ui,vi)} 

is  maximized  “over  one  zero”  p\  =  U\  +  jv ±.  Hence,  the  “Fano  1-bound” 
in  the  title  of  the  plot.  The  VSWR  is  computed  from  the  optimal  K  using 
Equation  12.  The  fractional  bandwidth  A /  determines  the  frequency  band 
[/min,  /max]  around  the  resonance  fr  =  843  MHz  as 

/min  =  /,(l-A//2), 

/max  =  /,(l  +  A//2). 

For  comparison  with  the  preceding  plots,  the  fractional  bandwidth  A /  = 
0.5932  approximates  the  interval  [600, 1100]  MHz. 


0.05  0.1  0.15  0.2  0.25  0.3  0.35  0.4  0.45  0.5  0.55 

fractional  bandwidth.;,/  =843  MHz 


Figure  19:  Microstrip  patch  antenna — Fano  1-bound. 
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Optimization  of  K  over  multiple  zeros  (u3,  v3,  u2,  v2l  u3,  v3,  . . .)  did 
not  change  the  plot  and  leads  to  the  natural  conjecture  that  Figure  19  is  the 
Fano  bound.  To  address  this  conjecture,  K(ui,vi)  is  plotted  in  Figure  20. 
The  plot  is  clipped  to  only  show  positive  K  values.  The  plot  shows  K(u3,vi) 
is  maximized  on  the  v\  axis. 


Figure  20:  Microstrip  patch  antenna — K  plot;  600-1100  MHz. 


Examination  of  the  functions  h0(u,v),  /ii(u, v),  h2(u, v),  and  h3(u,v) 
near  the  v  axis  is  listed  in  Table  5.  The  first  numerical  column  evaluates  the 
functions  on  the  v  axis  and  shows  all  functions  are  constant  for  u  =  0.  The 
last  column  lists  gradient  of  each  function  on  the  v  axis  with  v  >  0.  For 
clarity,  the  frequency  scaling  is  omitted  and  these  normalized  functions  are 
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denoted  by  the  capital  UH” : 


Hq{  u,  v) 
Hi(u,  v) 
H2(  u,  v) 


H3(  u,  v) 


Ur 


—3 A3  +  2 ^  '  um(um  —  3vm), 

-4? -2  E-ttt-. 

'  7/z  -U 
1 

-3^  +  2]T 


m(«m  -  3^) 


(«m  +  Vm)3 


This  column  shows  that  every  component  of  the  gradient  in  the  neighborhood 
of  the  v  axis  is  decreasing  in  the  u  directions.  Because  the  functions  are 
constant  on  the  v  axis,  this  decrease  forces  the  v  axis  to  be  a  local  maximum 
of  the  K  function.  Because  each  component  of  the  gradients  are  the  same,  the 
same  shape  of  K(u\,V\)  of  Figure  20  is  replicated  when  plotting  K(u,bfv). 
Consequently,  the  Fano  1-bounds  of  Figure  19  are  the  Fano  bounds. 


Table  5:  Fano  functions  near  the  v  axis. 


MO,  v) 

dHn(0,v)/dum 

ho 

0.2545 

-2 

hi 

0.2743 

— 6v2 

hr2 

0.2444 

—2v~2 

h:i 

0.3814 

— Qv~ ^ 
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8  Ideal  Dipole 


If  VSWR  is  plotted  as  a  function  of  frequency,  local  minima  identify  local 
resonances.  At  each  resonance,  the  local  sublevel  sets  are  intervals: 


[/i,/2]  =  {/eM:VSWR(/)  <5}. 


Consequently,  one  definition  of  bandwidth  A /  is  based  on  these  sublevel  sets 
of  the  VSWR  [21,  Eq.  5]: 


—  ■)  *  r  5 

Jr 

where  fr  is  a  resonant  frequency  in  the  frequency  band  [/i,  fo\. 

Hujanen,  Holmberg,  and  Sten  obtain  an  upper  bound  on  the  bandwidth 
A /  for  a  small  antenna  with  radius  a  and  wavenumber  k  [15,  Eq.  4]: 


A/< 


( ka )3 

1  + (ka)2 


S2  -  1 
25 


(18) 


provided  the  input  impedance  is  modeled  by  a  single  pole.  Figure  21  shows 
two  schematics  of  the  impedance  models  for  an  ideal  dipole. 


Figure  21:  Impedance  models  of  the  TMi  and  TEi  modes  of  the  ideal  dipole; 
C  =  ae o  =  a  x  8.85  (pF),  L  =  a/i o  =  a  x  1.26  (/xH),  R  =  \J p® /e$  =  377.3233 
(ohms). 
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Resonance  for  TMX  occurs  when  \zl{juS)\  is  minimum  [7,  page  81]: 


R(L  +  RCl/2(R2C  +  2  L)1/2)1/2 
Ur  ~  IM\R?>C2  +  2LR2  -  L 2)V2 

1  (1  +  31/2)1/2 
a  2 1/2v//i0e0 

=  ~(1  +  31/2)1/22-1/2. 
a 

In  terms  of  frequency, 

55.70 


Figure  22  plots  the  impedance  of  the  TMi  mode.  Figure  23  plots  the  associ¬ 
ated  VSWR.  Both  plots  show  that  bandwidth  is  a  local  definition.  The  pole 
at  /  =  0  and  the  flat  rolloff  at  the  high  frequencies  distort  the  symmetry  of 
/i  and  /2  with  respect  to  fr.  Figure  24  plots  the  dipole’s  reflectance  on  the 
Smith  chart. 


TMj  Impedance 


Figure  22:  Ideal  dipole — impedance  of  the  TMi  mode. 
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TMj  VSWR;  -0=50  ohms 


Figure  23:  Ideal  dipole — VSWR  of  the  TMi  mode. 


Figure  21  shows  that  the  shunt  inductor  is  a  short  circuit  at  p  =  0  and 
that  the  series  capacitor  is  an  open  circuit  at  p  =  0.  Theorem  4  identifies 
these  two  cases  as  transmission  zeros.  Because  the  degree  of  the  circuit  is  2, 
these  are  the  only  transmission  zeros.  The  reflectance  looking  in  from  the 
resistor  and  normalized  to  Rq  is 

,  ,  _  zx(p)  ~  Rq  _  zx(p)/R  -  1 
SX  P  zx (p)  +  Ro  zx (p)/R+  1  ’ 

using  the  notation  of  Equation  13  to  handle  the  resistor  R.  In  this  case,  the 
impedance  looking  in  from  the  resistor  is 

Zx -{p)  =  (C vL)~ 1  +  ((pC)-1  +z2(p))~1y1, 
where  Z2  (p)  is  the  impedance  looking  out  from  the  antenna  into  the  matching 
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Figure  24:  Ideal  dipole — reflectance. 


circuit.  Expanding  at  p  =  0: 


log  ( s  x  Go)  x) 


Ag  +  A"p  +  A"p6  +  A°p 4 

L  2L2(3R2C-L)  , 

-J7T  +  2— p - — - p 

+2j^pi+ob>_% 


provided  ^(p)  is  bounded  as  p  — >  0.  As  expected,  the  transmission  zeros  at 
p  =  0  pushes  the  load  22  (p)  into  the  higher  order  terms. 
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These  transmission  zeros  force  the  Fano  integral  constraints  [9,  Table  II] 

roo 

l  u-2  log(|sx(jo;)|_1)(i(n  =  -A\  -7 (19) 

/»oo 

Jq  ^“4  log(|sA-(jn;)|-1)^  =  -~A%  +  -  (20) 

where  the  pm’s  are  the  zeros  of  Sx(p )  that  lie  in  the  open  right  half  plane. 
Denote  the  frequency  band  as  D  :=  [camin,  wj.  Bound  the  performance  of 
the  reflectance  sx(p)  on  D  as 

min{log(|sx(jw)|-1)  :  uj  G  D}  =:  ttK/2. 

Consequently,  the  reflectance  Sx  (p)  is  bounded  on  the  entire  frequency  axis 
as 


^og(\sx(juj)\  X)  > 


— j  kK/2  camin  <  |oo|  <  oumax 
0  otherwise 


Inserting  this  lower  bound  into  the  Fano  integrals  produces  the  inequalities: 


7 iK 

- I 

2 

7 tK 


f^min  ^max ) 


< 


3 


3 


6 


^  ^min  ^max )  — 


07  2  log(|sA(j07)|  l)du  =  -A\  -  7 T^pJ, 
07'4  log(|sA (ju^-^du  =  ~A%  +  ^  ^p-3. 


Denote  the  zeros  of  Sxip )  as  prn  =  um  +  jvm.  Because  the  sums  are  real,  the 
Fano  integral  constraints  force  the  equalities: 


3*  [Y1 P™ 


,p 


-3 


E 

E 


Ur 


7/2  _j_  7,2  5 
m  1  m 

Um(u2m  ~  ?>V2m) 


(«m  +  Vl)3 

Combining  these  inequalities  with  the  Fano  integrals  produces 


K(u. 
K{  u. 


min  ^max ) 


<  A\-2 


Ur 


<  ^  77,2  _ L  77  2  ’ 

1  m 


3 


min  ^max )  — 


-3^  +  2jJ 


~ 

«  +  Um)3  ’ 
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The  goal  is  to  maximize  K  subject  to  the  inequalities  determined  by  the 
zeros  in  the  right  half  plane.  Define 


Vu,v)  :=  (c c, 


Um  -  3u>nVl 
(«m  +  VIY 


hi(n,  v)  :=  (tcmfn  -  uj^)  1  |-3AjJ  +  2  ^ 

The  Fano  bounds  are  computed  by  maximizing 

K  =  h0(  u,  v) 

subject  to  the  constraints  that  u  >  0,  v  >  0  and  that 

0  =  hi (u,  v)  -  h0( u,  v), 

where  the  “>”  is  used  to  simplify  the  discussion.  In  terms  of  Lagrange 
multipliers,  the  Lagrangian  can  be  written  as 

L{ u,  v)  :=  h0(u,v)  -  A{/ii(u,  v)  -  h0( u,  v)}. 

The  Lagrange  multiplier  theorem  identifies  local  extremum  as  stationary 
points  of  the  Lagrangian  or,  equivalently  [17,  page  243] 

(1  +  X)Vh0  =  AV/ii, 

provided  the  points  are  regular.  Because  each  vector  component  of  the  La¬ 
grangian  is  the  same  equation,  the  same  roots  apply  to  all  components  and 
the  resulting  four  roots  are  tabulated  in  Table  6. 


Table  6:  Extrema  of  the  Lagrangian  without  positivity  constraints. 


W  :  = 


Um 

Um 

v/3v/-A(l  + A)1F 

0 

-v/3v/-A(1  +  A  )W 

0 

0 

a/A(1  +  A  jw 

0 

—  a/A(1  +  A )W 

wn 


i2  +  wn 


F  1C min  HI; 


mm  "-  iriin 


W  n 


Wmin2A  +  1C 

max^  A  +  1C  max 2  +  W  max  ^min  ^  ^min  ^max  4“  ^mir4 
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What  is  fascinating  about  these  solutions  is  that  the  roots  are  parameterized 
by  the  frequency  band  and  the  multiplier — it  appears  that  the  load  does  not 
affect  these  roots.  However,  the  load  does  enter  in  the  roots  through  the 
multiplier  A  selected  to  enforce  the  constraint  equation.  The  positivity  of  u 
and  v  forces  two  cases  listed  in  Table  7. 

Table  7:  Positive  extrema  of  the  Lagrangian. 


u 

V 

A 

v/3v/-A(l  +  A)|W| 

0 

0 

y/X(l  +  X)\W\ 

-1  <  A  <  0 

A<— lorA>0 

However,  this  complicated  dance  through  the  inequalities  can  be  avoided 
and  more  insight  gained  by  considering  the  plot  of  K(ui,v\).  Figure  25 
corroborates  Table  7  by  showing  the  maximum  of  K(ui,Vi)  lie  on  both  axis 
with  the  maximum  on  the  V\  axis. 
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Because  the  components  of  the  gradients  V h0  and  V/ii  yield  the  same 
form,  Table  7  forces  the  maximum  on  the  v  axis  using  these  Fano  functions: 


AO 

V  0,v)  =  —  1 


^min  ^max 


^min  ^max 


-1  X  i?’ 


hi(Q,  v)  = 


34° 


^min  ^max 


_3  g 

^rnin  ^max 


X 


L2(3R2C  —  L) 

TR3 


Figure  26  compares  this  Fano  bound  for  the  VSWR  against  the  Fano  1-bound 
computed  by  numerically  optimizing  K(ui,vi).  Consequently,  the  analysis 
and  numerical  simulations  argue  that  this  plot  does  show  the  Fano  bounds. 


Ideal  Dipole:  Comparing  Fano- 1  bounds  and  analytic  solution 


fractional  bandwidth;./^  5 5 . 7 048  MHz;  a=  1 


Figure  26:  Ideal  dipole — Fano  bounds. 
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9  Automating  the  Fano  Bounds 


The  preceding  analysis  demonstrates  several  issues  must  be  addressed  to 
automate  the  Fano  bounds. 

Rational  Load :  Basic  to  the  Fano  analysis  is  that  the  load  is  a  rational 
function — not  sampled  impedances.  Mapping  sampled  data  to  a  bounded 
real  rational  function  is  a  substantial  problem  in  approximation  theory. 

Darlington's  Representation:  The  second  step  of  the  Fano  analysis  requires 
that  the  Darlington  representation  be  determined  from  the  load  to  compute 
the  scattering  matrix  Sx(p)  and  the  Fano  reflectance  sx(p)- 

Transmission  Zeros:  The  third  step  of  the  Fano  analysis  requires  the  trans¬ 
mission  zeros  of  Sxip)-  Computing  zeros  of  rational  functions  is  always 
tricky.  An  intermediate  approach  applicable  to  antenna  design  would  re¬ 
strict  the  analysis  to  transmission  zeros  at  p  =  0  and  p  —  oo. 

Fano  Functions:  The  fourth  step  is  to  generate  the  Fano  functions  h0( u,  v), 
hi(u,  v),  ...  associated  with  each  transmission  zero.  A  symbolic  manipulator 
can  compute  a  running  expansion  of  log(sx(p)_1)  and  halt  when  £2(p)  ap¬ 
pears  in  the  expansion.  A  program  could  generate  a  “Fano  file”  that  encodes 
each  of  these  functions  and  compute 


K  =  min{/r0(u,  v),  hi(u,  v), . .  .}. 


Maximizing  K:  The  fifth  step  of  the  Fano  analysis  is  to  maximize  K.  The 
preceding  step  has  written  the  Fano  functions  in  an  M-file.  Figure  20  shows 
the  usual  problems  with  finding  a  good  starting  point  and  crawling  to  a 
maximum  on  non-differentiable  surface. 

In  summary,  automating  the  Fano  bounds  appears  possible,  provided 

•  A  Darlington  representation  of  the  load  is  computable. 

•  High-quality  transmission  zeros  are  computable. 
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A  Fano-Bode  Integrals 


The  Fano-Bode  integrals  are  computed  in  [9,  Table  I].  This  appendix  re¬ 
derives  Fano’s  result  when  the  transmission  zeros  are  at  p  =  oo.  The  deriva¬ 
tion  reveals  Fano’s  adroit  choice  of  the  analytic  windows  that  peel  off  the 
reflectance  coefficients. 

Assume  the  Fano  reflectance  sx(p )  is  a  real,  rational  function  with  ex¬ 
pansion  at  p  =  oo  [9,  Eq.  9] 


A°° 

log  (sxip)-1)  =  j(3  +  —  + 

P 


p3 


+  ...+ 


/too  /loo 

/12V+1  /12V+2 


P 


2N+1 


P 


,2N+2 


where  the  Aj^’s  are  real  and  /3  is  either  0  or  —  n  depending  on  the  sign  of 
sxip)  at  infinity.  Assume  that  sx(p)  does  not  vanish  on  j’M.  Let  p±,  ..., 
Pm  denote  the  zeros  of  Sxip)  in  open  right  half  plane  C+.  The  Fano-Bode 
integral  is 


7 r 


w  n  log(|sA-(jw)|  )du  =  (-1)”- 


A°° 

/12n+l 


2  n  +  1 


„2n+l 

fm 


}• 


for  n  —  0, 1, . . . ,  N. 

Proof:  Define  the  Blaschke  function  [22]: 


B(p)  :  = 


n 

m=  1 


p-pm 
p  +  pm 


Because  the  zeros  are  in  the  open  right  half  plane  C+,  B(p)  is  analytic  on 
C+.  Because  the  zeros  are  either  real  or  come  in  conjugate  pairs,  B(p)  is 
real.  Define  e  =  signjlirrip^oc  s(p)]  and  set 

Sxip)  ■=  eBip^sxip). 


Observe  Sxip )  is  real,  analytic  on  C+,  and  has  no  finite  zeros  in  C+.  Consider 
integrating  over  the  contour  Cr  in  over  the  contour  of  Figure  27. 
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Figure  27:  Contour  Cr,  line  segment  Lr ,  half  circle  Hr. 
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Because 


MkvCMi  )  =  (-i) 


AO o 

\N  ^±2N+2 


UU 


2N+2 


+  • 


as  |u;|  — >  oo,  the  integral  converges  to  the  well-defined  integral: 


lim 


(juj)  nlog (sx(juj)  )du  =  (-1) 


r»oo 

\n  I  ,  ,2 n 


UJ  n\og(\sx(juj)\  )du 


for  n  —  0,  . .  .  N.  The  integral  on  the  half  circle  Hr  is  evaluated  using  the 
asymptotic  expansion: 


p2nlog(esx(p)  )dp 


'  Hr 


Ao o  400 

— +  ^  +  . 


400  400 
.  ^N+l  ^2^+2 

•  I  O  AT  1  1  I 


'Hr  L  P  P 

Consider  the  case  when  n  =  0: 
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2N+1 
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,2N+2 


dp. 


log  (esx(p)  )dp 


'Hr 


'Hr 


—  +  -4-  + 


p  p 

=  A °°  log  (p)  - 

=  -j*  AT. 

When  n  —  N: 


Aoo 
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Aoo  Aoo 
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Aoo 
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Gathering  the  results,  the  first  integral  on  the  right  of  Equation  21  decom¬ 
poses  over  the  contour  as 


pn\og(esx{p)  )dp 


'Cr 


=  (-1  y 


UJ  n  log(|sA-(iw)|  )du  -  jnA^n+1. 


(22) 


The  second  integral  on  the  right  side  of  Equation  21  decomposes  over  the 
contour  as 

[  p2n\og{B(p))dp  =  [  (. ju)2n\og(B(ju))jdu+  [  p2n  \og(B(p))dp. 

J  Cr  J-r  J  Hr 

Because  log (B{ju))  =  jaxg[B(ju)]  and  arg [g(ju)\  an  odd  function,  the  line 
integral  is  zero.  The  integral  on  the  half  circle  is  evaluated  from  the  expansion 
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Xlog 
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as  p  — >  oo.  Consider  the  case  when  n  =  0: 
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When  n  —  1, 


f  2 

/  p2l0g(S(p))dp  =  j-TT  V 

JHr  6  rn=  1 


Pm- 


When  n  =  2, 


M 


I Hr 


P 4  log(5(p))dp  =  J-7T  pjj, 

ra=l 


Gathering  these  results, 


r  9  M 

/  .  p2n  \og(B(p))dp  =  j  7T  P™+1. 

Jcr  m=l 


Substitute  Equations  23  and  22  into  Equation  21: 


f  oo  o 

.2 n  i - /l  „  /  •  ,\|-1\  - 


o  =  it-1)"  /  w  n  log(|s.Y(jw)|  )dw  -  jvrA^+1  +  j 


2n  +  1 
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££,+1 
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B  Scattering  Measurements  to  a  Darlington 
Representation 

Section  9  pointed  out  that  the  Fano  Method  requires  a  Darlington  represen¬ 
tation  of  the  load.  A  real-world  load  is  typically  measured  over  a  frequency 
band.  Extracting  a  credible  Darlington  representation  from  these  sampled  re¬ 
flectances  is  challenging.  The  state-space  representation  of  a  passive,  lumped 
iV-port  provides  a  physically  credible  model  automatically  disgorges  a  Dar¬ 
lington  representation.  Figure  28  motivates  the  discussion  by  approximating 
a  measured  antenna  reflectance  (blue  line)  using  a  state  space  model  (green 
line).  Although  the  finer  features  are  lost,  this  state-space  model  of  the  data 
used  only  two  inductors,  two  capacitors,  and  a  single  resistor  to  obtain  a 
worst-case  fit  over  2-30  MHz. 


sG=SCP  8FT  with  transformer;  2  </<  30  MHz 


Figure  28:  Surface  Current  Probe  approximation  with  degree  4. 
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Figure  29  illustrates  the  state-space  representation  of  a  passive,  lumped 
iV-port.  The  figure  shows  that  by  pulling  the  d  reactive  elements  and  the  r 
lossy  elements  into  the  augmented  load  S jJ  (p) ,  what  is  left  is  an  (N+d+r)- port 
consisting  of  the  wires,  transformers,  and  gyrators.  This  constant  (N  +  d  + 
r)-port  has  constant  scattering  matrix  Sa,  called  the  augmented  scattering 
matrix.  The  iV-port  is  obtained  by  looking  into  Ports  1,  2,  .  . . ,  N  of  the 
augmented  scattering  matrix  Sa  while  its  remaining  ports  are  terminated  in 
the  augmented  load  Sl{p).  That  is,  S(p)  is  the  image  of  the  augmented  load 
viewed  through  the  augmented  scattering  matrix. 


Figure  29:  State-space  representation  of  a  passive,  lumped  iV-port  containing 
d  reactive  elements  and  r  resistors. 

After  introducing  the  necessary  definitions,  Theorem  5  gives  the  precise 
statement  of  the  state-space  representation. 

Definition  1  Let  S  £  iR.BH°°(C+,CN><N)  and  be  rational  [26,  page  91]: 

•  The  normal  rank  of  S(p)  is 

r[S(p )]  :=  rank[/jv  -  S(-p)TS(p)]. 
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•  The  skew- symmetric  rank  of  S(p)  is 

g[S(p )]  :=  ^rank^p)  -  S(p)T}. 

•  Let  the  K  distinct  poles  of  S(p)  be  denoted  as  pk.  Let  deg(S(p);  pk) 
denote  the  largest  order  to  which  p  =  pk  appears  in  any  minor  of  S(p). 
The  Smith- McMillan  degree  of  S(p)  is 

K 

degSM[S(p)]  :=  J]deg (S(p)-pk). 

k= 1 


Theorem  5  (State-Space  [26,  pages  90—93])  Every  lumped,  passive,  ca¬ 
sual,  time-invariant  N -port  admits  a  scattering  matrix  S{p)  and  conversely. 
If  S  (p)  has  Smith- McMillan  degree  d  and  normal  rank  r,  then  S(p)  admits 
the  following  state-space  representation: 

S(p)  =  E(Sa,SL]P)  :=  Sa,n+Sa,l2SL(p)(Id+r  ~  S^Sl)'1  Sa,21, 


where  the  augmented  load  is 


Sl(p) 


q  x  INl  0  0 

0  —q  x  I  nc  0 

0  0  0  x  Ir 


5 


and  Nl  +  Nc  =  d.  The  augmented  scattering  matrix  is 


Sa 


Sa  ,11  Sa  ,12  N 
Sa, 21  Sa, 22  _  d  +  r 

N  d  +  r 


is  a  constant,  real,  orthogonal  matrix.  If  the  N -port  is  lossless,  then  r  —  0. 
If  the  N -port  is  reciprocal,  then  Sa  is  symmetric:  Sff  =  Sa. 


Although  a  scattering  matrix  S(p)  has  lots  of  representations,  the  fol¬ 
lowing  list  shows  that  the  state-space  representation  extracts  a  lot  of  the 
iV-port’s  structure  [26,  page  91],  [19].  For  brevity,  “in  any  passive  IV-port” 
means  the  A^-port  of  Theorem  5. 
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Resistors.  The  number  of  resistors  Nr{  in  any  passive  Af-port  exceeds  the 
normal  rank: 

Na>r[S(p)]. 

Moreover,  of  all  the  AT-ports  that  have  S (p)  as  their  scattering  matrix, 
there  exists  at  least  one  having  exactly  r[S(p)]  resistors. 

Reactive.  The  number  of  inductors  Nl  and  capacitors  Nq  in  any  passive 
Af-port  exceeds  the  Smith-McMillan  degree: 

Nl+Nc>  degSM[S(p)]. 

Moreover,  of  all  the  AT-ports  that  have  S(p)  as  their  scattering  matrix, 
there  exists  at  least  one  AAport  having  exactly  this  many  inductors  and 
capacitors. 

Gyrators.  The  number  of  gyrators  in  any  passive  Af-port  exceeds  the  skew- 
symmetric  rank: 

Ng  >  g[S(p )]. 

Wohlers  [26,  page  92]  remarks  that  only  in  the  case  that  g[S(p)\  =  0  or 
that  S(p)  is  symmetric,  can  an  Wport  can  be  found  that  is  reciprocal 
or  contains  no  gyrators. 


Theorem  5  provides  a  numerically  efficient  parameterization  of  the  scat¬ 
tering  matrices  corresponding  to  the  lumped,  passive  ./V-ports  independent  of 
any  circuit  topology.  Denote  the  class  of  lumped,  lossless  Wports  of  degree 
d  as 

U*(N,d)  :=  {S  €  U+(N)  :  degsM[S(p)]  <  d}. 

Next,  introduce  the  group  of  M  x  M  orthogonal  matrices: 

0[M\  ■=  {S  G  MMxM  :  STS  =  IM}. 

Theorem  5  determines  that  U+  (N}  d )  is  the  orbit  of  a  fixed  reactive  load 
under  the  action  of  the  orthogonal  0[N  +  d\. 


Theorem  6  Any  augmented  load , 


has  orbit 


Sl(p ) 


p  —  1 
p+  1 


INl  0 
0  —Inc 


(Nl  +  Nq  =  d ) 


U+(N,d)  =T(0[N  +  dfiSL). 
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Sweeping  over  the  orthogonal  group  0[N  +  d]  generates 
all  lumped,  lossless,  IV-ports  of  degree  d. 


It  is  known  that  U+(N,  d )  is  compact,  so  best  approximations  from  U+(N,  d) 
do  exist.  Unfortunately,  little  is  known  regarding  the  characterization  of 
best  approximations  from  U+(N,  d)  and  the  rate  of  convergence  as  d  — »  oo. 
Therefore,  the  remainder  of  this  section  simply  shows  the  quality  of  numer¬ 
ical  approximations  for  several  real-world  antennas.  The  goal  is  to  obtain 
a  Darlington’s  representation  for  a  given  reflectance  sG.  The  first  step  is  to 
compute  the  best  approximation  from  the  orbit  of  a  1-ohm  resistor  under  the 
action  of  the  lossless  2-ports: 

min{||sG  -  Si||n,oo  :  «i  e  U(N]  Nl,  Nc,  Nr)}, 

where  the  notation  returns  to  the  decomposition  of  Theorem  5.  The  number 
of  ports  is  denoted  by  N.  Because  the  reflectance  sq  is  a  scalar,  N  —  1.  The 
number  of  inductors  and  capacitors  is  denoted  by  NL  and  Nc,  respectively. 
The  number  of  resistors  is  denoted  by  Nr.  Because  a  Darlington’s  represen¬ 
tation  is  being  computed,  Nr  =  1.  The  Darlington’s  representation  takes 
the  reflectance  si  at  Port  1  with  the  understanding  that  the  resistor  is  lo¬ 
cated  at  Port  2.  The  plots  reports  the  reflectance  as  si  e  17(1,  Nr ,  Nc,  Nr). 
The  definition  of  “degree”  is  slightly  expanded  as  deg[si]  =  [Nr  Nc  1]  to 
make  explicit  that  a  the  single  resistor  is  used.  The  caption  to  each  plot  still 
reports  the  degree  d  as  the  number  of  reactive  elements. 
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Figures  30  and  31  show  the  quality  of  this  approximation  to  the  re¬ 
flectance  of  the  AOC21B  measurements  [23].  The  reflectance  is  plotted  in 
blue  and  shows  two  large  resonances  and  dozens  of  small  resonances.  The 
approximations  are  the  green  lines  and  a  visually  appealing  fit  is  obtained 
with  a  16-degree  2-port. 


sG= aoc21  b;  2  </<  8  MHz 


Figure  30:  AOC21B  and  approximation  with  degree  8. 
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sQ='c U0C21  b;  2  </<  8  MHz 


Figure  31:  A0C21B  and  approximation  with  degree  16. 
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The  CNV77Z  measurements  (blue  line)  are  less  complex  than  the  preced¬ 
ing  AOC21B  measurements.  Relatively  good  approximations  are  obtained 
for  small  degree  d  as  shown  in  Figures  32  and  33.  Nevertheless,  Figure  34 
shows  that  the  small  resonance  is  not  approximated  until  degree  16  is  reached. 


sr  =CVN77Z;  2  </<  9  MHz 

G  ’  J 


Figure  32:  CVN77Z  approximation  with  degree  4. 
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sg=CVN77Z,  2  </<  9  MHz 


Figure  33:  CVN77Z  approximation  with  degree  10. 
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sg=CVN77Z,  2  </<  9  MHz 


Figure  34:  CVN77Z  approximation  with  degree  16. 
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The  Surface  Current  Probe  started  this  discussion.  Figure  28  shows  the 
fit  when  the  antenna  has  been  “pre-processed”  by  a  transformer.  The  trans¬ 
former  pnlls  the  reflectance  further  into  the  Smith  chart  a  fine  structure 
appears.  Comparison  with  Figure  35  is  interesting  because  both  plots  have 
the  same  degree  and  the  measurements  differ  only  by  a  transformer.  If  the 
transformer  was  ideal,  the  errors  ||sg  —  si||oo  should  be  identical  because  the 
lossless  2-ports  sweep  over  all  transformers — all  ideal  transformers.  However, 
the  real-world  transformer  induces  resonances  that  prevent  a  tight  fit. 


v,  =SCP  8FT;  no  transformer;  2  </<  9  MHz 


Figure  35:  Surface  Current  Probe  (no  transformer)  approximation  with  de¬ 
gree  4. 
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